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Abstract 

We measured the neutron time-of-flight spectrum in the 4 He(stopped-fr _ , n) reaction 
by stopping negative kaons in a superfluid helium target. A clear enhancement 
was observed in the neutron spectrum, which indicates the formation of a strange 
tribaryon of charge +1 with a mass and width of M s + = 3141 ± 3 (stat.) ^ 4 
(sys.) MeV/c 2 and T s + < 23 MeV/c 2 . This state, denoted as S+(3140), is about 25 
MeV/c 2 higher than the previously observed S°(3115) T = 1. 
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1 Introduction 



Very recently, we have discovered a strange tribaryon, S°(3115), via the 

(K~ 4 He) 

atomic 

S u (3115)+p (1) 

reaction [1]. The observed state has isospin T — 1, charge Z = and mass 
M s o ~ 3117 MeV/c 2 . 

Where is the isospin T = state? It is also important to search for an isospin 
partner of the observed state, S°(3115) T — 1. The S°(3115) discovery was 
originally motivated by the theoretical prediction of a deeply-bound kaonic 
state by Akaishi and Yamazaki at M — 3194 MeV/c 2 with T = and Z — 1 
[2], which can only be studied via the 

(K~ 4 He) atomic -> S+ + n (2) 

reaction. In this reaction, both isospin T = and 1 states can be populated, 
while the reaction (1) is limited to forming T — 1. We describe the study of 
the neutron spectrum obtained from reaction (2) to provide further important 
information on the strange tribaryon. 



2 The experiment 

Let us briefly summarize the experiment while focusing on a neutron mea- 
surement. A more detailed description is given in Refs. [1,3]. 

To obtain the neutron spectrum, we adopted the time-of-flight (TOF) method, 
as shown in Fig. 1. For this purpose, we need to define the kaon reaction 
point by the vertex between the trajectory of an incoming kaon and that of 
an outgoing charged particle detected by two layered trigger counters (TC), 
namely TC t hi n and TC t hi C k- The resolution of the reaction points/vertex was ~ 
5 mm (a), which is predominantly determined by kaon multiple scattering in 
the target helium. The kaon reaction at-rest in the target volume was selected 
(purity > 95%) by using the correlation between the pulse height of the final 
beam counter (TO) and the range in the helium target obtained by the vertex 
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Fig. 1. Schematic figure of the experimental setup. The notations are: TO, beam 
timing counter; BDC, beam-line drift chamber (16 layers); VDC, vertex drift cham- 
ber (12 layers); TC t hi n and TC t hick> thin (0.6 cm) and thick (3 cm) trigger counter; 
NCV, neutron-counter charged-particle veto; NC, neutron-counter wall. To simplify, 
several counters are omitted. 



position. The reaction (neutron TOF start) timing was calculated from the 
TO timing and the kaon range in the target, and the neutron was detected 
by two sets of plastic scintillation counter arrays (NC). An absolute time 
calibration of the TOF was performed by electromagnetic showers produced 
in iron-plate converters, placed between NC veto counters and NC, itself. The 
overall TOF resolution was measured by this shower event, and was found to be 
A(l//5) = 0.04(it). The absolute momentum was checked using mono-energetic 
neutrons from S + decay at rest (S + — > nn + ) by back-to-back coincidence with 
the threshold at 3 MeV electron equivalent (MeVee) applied to NC. 



In a following analysis, we set a rather high energy threshold to the NC system 
(10 MeVee) at the analysis stage, so as to reduce the accidental neutron back- 
ground. At this threshold, an NC counter of 5 cm thickness become sensitive 
at ~ 150 MeV/c, and saturate at ~ 340 MeV/c. Beyond this, the efficiency is 
constant at ~ 20%. 
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3 Analysis procedure 



As shown in Fig. 2 (upper left), the semi-inclusive neutron spectrum does not 
yield any distinct peak, in contrast to the proton spectrum [1]. 
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Fig. 2. Neutron momentum spectra. The upper-left panel shows the semi-inclusive 
spectrum without any event selection, whereas the rest of the figures are subsets. 
The event selection of each subset is given in Fig. 3. A peak- like structure exists in 
a "fast 7r" triggered event set at around ~ 470 MeV/c; a close-up view of the doted 
region is given as an inset. 

To study the neutron spectrum in more detail, we focused on pions that came 
from hyperon decay, focused on pions that came from If these pions come from 
hyperon decay, the mother hyperon could be qualitatively identified by the 
pion momentum. The momentum of the pion from A distributes dominantly 
in the lower momentum side (mostly below 100 MeV/c) compared to that 
from a charged S (centered at about 190 MeV/c) [4]. 

For a further study, we proceeded to an analysis of the hyperon motion based 
on the vertex inconsistency. If the hyperon moves a certain length, the two 
trajectories between the incoming kaon and the outgoing charged pion could be 
separated by the motion before the decay. On the other hand, if the triggered 



4 



pion come from the quasi- free hyperon production (QF : K N — > ir^Y), the 
vertex should coincide within the experimental resolution. 

If the main decay channel of S + is similar to that of S°(3115), then one can 
assume that the EiVTV is the major decay channel. Thus one expects 

S+ _> Z±NN and E ± -> ^N, (3) 

which is an ideal model to apply the analysis procedure described above, 
because S + is expected to be boosted at the formation stage. 



4 Neutron spectra classified by TC 

Let us define the event window using the TC information. Fig. 3 shows the 
scatter plot of dE/dx on two layers of the TC counters. To select high- 
momentum pions (p-x 125 MeV/c), which mainly come from S 1 * 1 decay, 
we proceeded as follows. From the simulation, we know that the pions, which 
penetrate both layers (p w ^ 100 MeV/c), should distribute mostly inside the 
dashed region. We evaluated the distribution center by using an average of the 
data within this region. The central track obtained is shown in the figure as a 
solid line, and the evaluation of the momentum along the track led us to select 
the "fast 7r" triggered event window (hatched area in the dashed region). This 
window is well apart from the "proton" window shown in the figure. 

The "slow 7r" triggered event window is defined so as to select pions that 
stop in the TC counters (p w ^ 90 MeV/c). Apart from the region shown in 
the figure, we also included a very slow pion which had already stopped in 
TCthin, because these events can be easily selected by the pulse height and 
TDC information. 

Fig. 2 shows the neutron momentum spectra for "proton" (lower left), "slow 
7r" (upper right) and "fast 7r" (lower right) triggered events tag. There exists 
a common structure up to ~ 400 MeV/c. This low-momentum component 
mainly comes from the 7r~ absorption reaction (occurring in or around the 
target), ir~NN — ■> NN [5], and cascade reactions followed by QF-hyperon 
formation [3]. 

The "slow 7r" triggered neutron spectrum has a "plateau" structure beyond ~ 
450 MeV/c. This comes from the kaon non-mesonic absorption and its chain, 

K~NN -> An AN and A -»■ n'p, (4) 

where uan is the neutron from the non-mesonic kaon absorption reaction. Note 
that pions from E decays do not contribute when we select low-momentum 
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Fig. 3. Scatter plot of dE/dx on TC t hick (horizontal-) and TC t hin (vertical-axis) 
counters. Definition of "proton", "slow 7r" and "fast 7r" triggered events are given 
in the figure. 

pions [4]. The spectrum shape of is dominantly determined by three- 
body phase space and the nuclear form factor, and is very broad, as shown 
in the figure. It has a maximum momentum of ~ 700 MeV/c, and forms a 
"flat-plateau" below ~ 600 MeV/c. 



Similar to the "slow 7r" triggered events, the neutron spectrum of "fast 7r" 
triggered events is expected to have a high-momentum component due to 

K~NN -» E^nsnipi p) and E ± -> 7r ± n SD , (5) 

where n^N is the neutron at the primary reaction and n-^D is that from 
decay. The n^N should have a nature very similar to the tt-an in reaction (4), 
but shifted by ~ 100 MeV/c to the lower side. Therefore, a "plateau" below ~ 
500 MeV/c is expected. Actually, the "proton" triggered event is not selective 
to n\ N , n SN or n^D, as indicated in Fig. 2 (lower left). 

There exists a peak-like structure at ~ 470 MeV/c in the "fast 7r" triggered 
event (Fig. 2 (lower right)), at the momentum where we expect a "plateau" 
due to n SN and ti-^d- Is this peak-like structure due to the formation of the 
intermediate state S + ? As shown in the figure, the answer largely depends on 
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what the background state is. and where the background level is around this 
region. 

If this peak-like structure indicates the formation of S + , then the major decay 
will be EiViV, because the structure appears only in the "fast 7r" triggered 
spectrum. To make the existence clear, we need to study the spectrum in 
more detail by analyzing the hyperon motion. 



5 Neutron spectra selected by the hyperon motion 

The most simple way to detect the hyperon motion is by a vertex inconsis- 
tency between the two tracks of an incoming kaon and an outgoing trigger 
particle (so called distance-of-closest-approach). In the present case, however, 
the hyperon and the neutron motions are asymmetric in many reactions (see 
Fig. 4). For a further analysis, let us introduce a scalar product, v CA -v n , where 
vca is the vector from the kaon trajectory to that of the charged pion at the 
minimum distance, and v n is the normalized vector of the detected neutron. 

Fig. 5 shows the v C a • v n distribution of "fast 7r" triggered events having NC- 
detected neutron momenta p n > 400 MeV/c. For a comparison, the decom- 
posed spectra obtained by a simulation are shown. As shown in the simulation, 
n ZN + ^an and n^D roughly distribute symmetrically on opposite sides. 

If the peak-like structure obtained in Fig. 2 (lower right) is due to the forma- 
tion of an intermediate state, S + , then it should be boosted opposite to the 
neutron with the same momentum, as shown in Fig. 4 (top). It should then 
decay dominantly to ENN rather than ANN, as is the case for S°, since we 
observe the structure only for a "fast 7r" triggered event. What is the centroid 
of the momentum of the decay E ±r ! We calculated the momentum, and found 
that it is about 450 MeV/c in the backward direction, while it is about 80 
MeV/c in forward direction. Thus, the vca ■ v n distribution for such events 
should be asymmetric and dominantly located on the negative side. This is 
the reason for the shape of the distribution of the simulated signal given in 
Fig. 5. 

Unfortunately, there is no clear criterion about how to separate the signal from 
the background. Simple "forward" and "backward" comparisons may not be 
valid, because the yield and the spectral shape can also be weakly asymmetric 
between the "forward" (tied) and "backward" (n^N + «an) directions. 

Therefore, we examined how to compare the event sets between higher S/N 
(HSN : indicated with arrow) and lower S/N (LSN : shallow hatched) regions, 
whose windows are defined as shown in the figure. HSN is defined to be where 
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Fig. 4. Example of the event topologies. The S + formation topology is shown at the 
top. If the Q-value of S + formation is large, then the S + is boosted to the opposite 
direction of the signal neutron, which results in the vca • v n being distributed on the 
negative side. The background reaction (5) is shown at the bottom. Although the 
two topologies are quite similar, the scalar product, vca ■ v„, is negative (backward) 
if nsN is detected by the neutron counter, while it is positive (forward) for tied- 



we expect a better signal-to-noise ratio than in the other regions. The reason 
why we sandwiched HSN in between the two LSN regions (inner and outer 
side), is to obtain an interpolation of the shape of the spectrum. It should 
be noted that both the spectrum shape and the statistics of the inner and 
outer LSN regions are different. If vca ■ v n is close to zero, then the pion 
dominantly originates from the primary QF reaction, namely K~N — > it^Y . 
If it is much apart from zero, the pion is dominantly from higher momentum 
hyperon decay, either from reaction (4) or (5), which results in larger yields 
of the higher momentum component in the neutron spectrum. 

It should also be noted that the LSN data are absolutely not free from the 
signal. In the case of the neutron spectra, the ratio of the signal yields between 
HSN and LSN is expected to be about ~ 2 : 1 from the simulation, if we assume 
that ^NN is the major decay branch of S + . 



8 



t- "fast Jt" triggered A 

p > 400 MeV/c 



10" 



simulated signal 
_ S + ->£" pp 

simulated f\-£ N + f[^ N 



§ 10 

o 

o 



10 




data 

simulated neutron from 
QF-ic triggered event 



simulated PI^d 



\ lower S/N region " JL. 

\{inner side) 



50 
V CA ■ V n (mm) 



Fig. 5. vca • v n distribution. The dots with error bars show the data of the "fast 
7r" triggered events associated with a high- momentum neutron {p n > 400 MeV/c). 
The decomposition into key components, given by the simulation, is also shown. To 
make the signal formation component visible, the S + is simulated with a formation 
probability of 1% per stopped K~ 0, and decays purely into T,~pp, which is the most 
efficient branch to detect in this vca ■ v n selection. 

6 Discussion and Conclusion 



Fig. 6 shows a comparison of the two spectra. It exhibits a clear enhancement 
at around ~ 470 MeV/c in the spectrum of HSN, whereas it is smooth in LSN. 
Because the windows for HSN and LSN are rather arbitrarily defined, we also 
examined the stability of the enhancement by changing the cut region, and 
verified that it was nearly unchanged. The same spectra-shape comparison, 
whose windows were chosen to be symmetric to the present one in the vca • v„ 
positive side, was applied to check the validity of the comparison. We found 
no clear excess, as was expected. 

The missing-mass spectra both on HSN and LSN, obtained from the neutron 
data, are given in Fig. 7 (bottom). It shows that an enhancement exists at 
around the mass region from 3110 to 3160 MeV/c 2 . As a third and final ex- 
amination, we applied exactly the same procedure to obtain the mass spectra 
from proton data, as shown in Fig. 7 (top). We fitted the data of the HSN, 
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Fig. 6. Comparison of the momentum spectra between HSN (closed circle) and 
LSN (open circle). The LSN data were scaled so as to have the total yield be the 
same as that of the HSN data, including the lower momentum region. To avoid the 
error bars from overlapping, the center of the bin is shifted by half of the bin size. 
The peak-like structure, observed in the "fast 7r" triggered neutron spectrum (Fig. 2 
(lower right)), appears as a clear excess between HSN and LSN. 



assuming a smooth background for both the proton and neutron spectra. The 
fit for the proton HSN spectrum (background shape of a third-order polyno- 
mial was used) gave a consistent result with that obtained in a previous paper 

[1]- 



The excess in the neutron spectrum seems to be rather broad compared to that 
in the proton spectrum. It is natural to expect the isospin partner of S°(3115) 
in this energy region; hence, we fitted the HSN neutron spectrum with two 
Gaussian functions together with a smooth background (a single exponential), 
as shown in the figure, while fixing the mass and width of one Gaussian to 
be the same as that of S°(3115). In the fit, data of LSN were not utilized 
to constrain the background shape/yield. Because the signal is expected to 
appear in both HSN and LSN at a ratio of 2 : 1, the background is expected 
to be located slightly below the LSN data. The result, however, is close to the 
data point of LSN, so that the fitted background slightly overestimated the 
yield in the region of interest. 
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Fig. 7. Missing-mass spectra of the 4 He( s t p P ed-^ _ )P) reaction (top) and the 
4 He (stopped-^ - ; n ) reaction (bottom). Both proton and neutron HSN spectra were 
fitted assuming a smooth background. In the case of the HSN neutron spectrum 
(bottom), we couldn't achieved a good fit without assuming two Gaussian functions 
in the region of interest. 
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It is obvious that the statistical significance of the isospin partner of S°(3115) 
is not enough to conclude its existence from the data itself. From the isospin 
rule, one can expect the formation ratio of the T = 1 states between the 
neutron and the proton emission processes to be 1 : 2. The lower neutron 
detection efficiency makes the ratio to be about 1 : 10. If we assume that the 
decay branching ratio to the charged S is the same from the two states, we 
expect the signal yield in the neutron spectrum to be about 1/10 of that in 
the proton spectrum, which is consistent even with the disappearance in the 
neutron spectrum. 

On the other hand, concerning the main signal component, the fit gave an 
energy of 3141 ± 3 MeV/c 2 , width of 17 ± 5 MeV/c 2 and a yield of 120 ± 
32 counts. By folding the systematic error and the experimental resolution, 
we obtained a mass and width of M s + = 3141 ± 3 (stat.) ±\ (sys.) MeV/c 2 
and i~s+ < 23 MeV/c 2 (95% CL), respectively. The significance of this excess, 
evaluated from the area of the Gaussian signal and its error, is 3.7 a. The 
observed excess indicates the formation of another strange tribaryon (denoted 
as S + (3140) hereafter), which lies about 25 MeV above S°(3115). Because of a 
lack of information about its decay branching ratio, we can not derive a reliable 
formation probability of S + from the fitted yield. Based on simulations, the 
formation probability is below a few % per stopped K~; otherwise, the peak 
should already be observed in the semi-inclusive neutron spectrum. 

Since the mass of the signal is much different from S°(3115), and also the yield 
is larger than the expected one, it is unlikely to be its isospin partner. This 
state could be assigned to the predicted T = state, although the observed 
mass is much lower than the theoretical one (3194 MeV). 

In summary, we have found an experimental indication of another strange 
tribaryon S + (3140) in the 4 He(stopped--ft'", n) reaction. Its isospin is likely 
to be T = 0, and the major decay mode is S + — > H^NN. The present find- 
ing provides additional information about the strange tribaryon system. The 
most important finding is that both S°(3115) and S + (3140) have much smaller 
masses (namely, much larger K binding energies) than the theoretical values. 
If these strange tribaryons, S°(3115) and S + (3140), are indeed deeply bound 
kaonic states, then the huge total binding energies (~ 200 MeV) naturally 
imply the formation of a cold and dense nuclear system [2,6]. These findings 
provide important information about the KN interaction, and may provide 
rich information to understand the dynamics of a dense system, such as a 
strange star [7]. In such a system, a restoration of chiral symmetry can be 
realized, and quarks may be deconfined [8], which are more or less unknown 
as of today. 
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